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Direct and Related Far-Field Inverse 
Scattering Problems for Spherical 
Electromagnetic Waves in Chiral Media 

N.M. Berketis * 

Independent Researcher, Athens, Greece 

Abstract 

This paper studies the direct and inverse scattering 
problem when the incident electromagnetic field is a time 
harmonic point - generated wave in a chiral media and the 
scatterer is a perfectly conducting sphere. The exact 
Green's function and the electric far-field patterns of the 
scattering problem are constructed. For a small sphere, a 
closed-form approximation of the scattered wave field at the 
source of the incident spherical wave is obtained. Also 
treats the same inverse problem using far-field results via 
the leading order term in the low-frequency asymptotic ex- 
pansion of the scattering cross-section. 



Introduction 

In a homogeneous isotropic 
chiral media the electromag- 
netic fields are composed of 
left - circularly polarized 
( LCP ) and right - circularly 
polarized (RCP) components, 
which have different wave 
numbers and independent di- 
rections of propagation. 

The LCP and RCP components 
are assumed to be spherical 
Beltrami fields since in 
practice such wave fields are 
more readily realized. 

In this work, the author 
has studied the electromag- 
netic waves in chiral media 
produced by a point source in 
the vicinity of the scat- 



terer. In particular, in [1], 
[2], reciprocity, optical and 
general scattering theorems 
for stimulation of point- 
source asymmetric media have 
been demonstrated. This paper 
studies the inverse problem 
of far field [2] . Specifi- 
cally, we measure the scat- 
tering cross-section for a 
five-point source area. 

In the second Section, 
considering Bohren decomposi- 
tion into suitable Beltrami 
fields, we formulate the di- 
rect scattering problem of a 
spherical electromagnetic wa- 
ve by a perfectly conducting 
obstacle. This problem is 
well posed, the existence and 
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uniqueness has been proved in 
[1], [3]. 

In the third Section, af- 
ter expanding the incident 
field in terms of spherical 
wave functions, we obtain the 
exact solution of the scat- 
tering problem as well as an 
expansion for the electric 
far-field pattern [2]. 

Finally in the fourth Sec- 
tion, we consider either LCP 
or RCP incidence and we ob- 
tain an approximation of the 
scattering cross-section. For 
the far-field experiments, we 
measure the scattering cross- 
section for various point 
source locations [2] . 



vector r 0 with respect to an 
origin 0 in the vicinity of 
the scatterer. In order to 
define spherical electromag- 
netic fields E r , H, , we make 
use of the Bohren decomposi- 
tion into Beltrami fields 
Q t , r , and , as follows 

E r e = Qt,r s + Q»,r a 

1 ( 1 ) 

H, ~ (Ql-,r 0 " Q/?,r 0 ) 

Where q = (p/s) 1/2 is the in- 
trinsic impedance of the 
chiral medium. The Beltrami 
fields satisfy the equations 
[4], [5], 



Statement of the problem 

Our goal is to study the 
direct and inverse scattering 
problems when the incident 
electromagnetic field is a 
time harmonic point - gener- 
ated wave in a chiral medium 
and the scatterer is a per- 
fectly conducting sphere of 
radius a centered at the 
origin. The exterior space 
(r = |r| > a) is an infinite ho- 
mogeneous isotropic chiral 
medium with chirality measure 
(3, electric permittivity s 
and magnetic permeability p. 

We consider a time har- 
monic spherical electromag- 
netic wave due to a point 
source at P 0 with position 



Vx Ql.,r 0 ~ K<-Qf-,r 0 

^ X Q/?,r 0 = — V'rQi.,!-,) 



where y L andy R are wave num- 
bers given by, 



Yl l - kp ' Yr l + fcp 



with k = <y(sp) 1/2 , w being the 
angular frequency. The indi- 
ces L and R denote the LCP 
and RCP fields respectively. 
The spherical incident Bel- 
trami fields with suitable 
normalization have the form 
[ 1 ], [ 2 ], 
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Qi n r :(r|pJ = -(I + — VV + — Vxi) ( 

z y l y \ 



K (Yl |r ~r B |) 
K (y L r b) 



2nr g 1Yl ~ rQ ~ ~ 

£ [y L G fs (r, r 0 ) + VxG fs (r, r 0 ) ] • p L 



Q^(r|p R ) = -(i + - 5 -W- — Vxi) ( 
^ Yr Yr 



K (Yr |r -r B |) 
h 0 (y R r B ) 



2tt r e 1>v " e 

£ [y R G fs (r, r 0 ) - VxG fs (r, r 0 ) ] • p R 



(4) 



(5) 



where b 0 (x) = b 0 (x) = e lx /(ix) is 
the zeroth-order spherical 
Hankel function of first 

kind, I = xx + ijjijj + zz is the 
identity dyadic, r 0 =|r o | and 

G fs (r, r B ) is the free space 
dyadic Green function [4] . 
The constant unit vectors 
p L and p R satisfy the rela- 
tions 

V Pl= *V P>?= o 

r 0 xp L = ip, (6) 

r 0 xp R = -ip R 

We note that when r 0 -»oo, the 
incident electric field 

(T | p t , p w ) = 

( 7 ) 

= Qt"r" ( r I PJ + Qr™ (r I Pr) 



reduces to plane electric 
wave with direction of propa- 
gation -r 0 and polarizations 
P L , p R , since 

limQf-(r|pJ = e 1 ^ r p,= 

( 8 ) 

= Qf c (r;-f 0 ,p L ) 

limQ™(r |p R ) = e-^ r p R = 

r ^° ° (9) 

= Qf c (r; -r B , p R ) 

We consider E“ c is incident 

r 0 

upon a perfectly conducting 
sphere of radius a . Then, we 
want to calculate the scat- 
tered electric field E* c , 

which is the unique solution 
of the following exterior 
boundary value problem 
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( 10 ) 
( 11 ) 

The Silver-Muller radiation condition is modified as follows 
r x VxE* c (r) - Py 2 r xE* c (r) + E* c (r) = o(y ) , r ® (12) 



uniformly in all directions 
reS 2 , where S 2 is the unit 

sphere in M 3 , n is the out- 
ward normal unit vector on 

the scatterer and Y 2 =Y L Y R - 
The scattered electric field 
will be depended on the po- 
larizations p L , p R and will 
have the decomposition 

E r 0 C ( r I Pl , Pr) = Q.1% (r I Pz., Pr) 

. . (13) 

+ Q"r 6 (r I P„ Pr) 

where Q“ r< (r | p L , p R ) and 



The functions g L and g Rr 
are the LCP and RCP far-field 
patterns respectively [4], 
[ 6 ], 

If either a LCP or a RCP 
spherical electric wave 
Ef" c (r|pj, A = L, R , is inci- 
dent upon the scatterer, then 
the scattering cross-section, 
is given by [6], 



f lA |g,,r 0 (r|pj| 2 + 

s 2 Yl 

+ A |g R ,r 0 (r I pjf] ds{r) (16) 
Yr 



Q“r e (r I P„ Pr) are the corre- 
sponding scattered Beltrami 
fields which have the follow- 
ing behavior, when r oo 

^{r\p Ll p R ) = h.{ YL ry 

(14) 

•g L ,r 0 (r |p L ,p R ) + 0(^) 

Q"r e ( r IP^P R ) = h B (Y R r)- 

. . . (I 5 ) 

g R ,r 0 (r|P^PR) + 0(^) 



Exact Green's function 

We take spherical coordi- 
nates ( r,8,cp ) where 9 e [0, n] 
and (p e [0, 2tt ) , with the ori- 
gin at the center of the 
spherical scatterer, so that 
the point source is at 

r = r B , 0 = 0 . 



Thus, r 0 



= r 0 z , P,= 



72 »-**> 



and p R = —j= (x + iijj) , where x, ip 
v 2 



and z are unit vectors in 



FUNKTECHNIKPLUS # JOURNAL 



V2-10 



ISSUE 2 - YEAR 1 




DIRECT AND RELATED FAR -FIELD INVERSE SCATTERING PROBLEMS 



the x, l/j and z directions, 
respectively. Using spherical 
vector wave functions, [4], 
[5], [7], and taking into ac- 
count (4), and (5), we obtain 
[ 2 ], 

oo 

Qt,™ ( r I Pi.) = X { L oln(/l.r) + 

->=1 (17) 

+ iLW„(/ t r)} 



Q^(r|p R ) = Xe:{R-(/ R r)- 



-iRW(y R r)} 



where, for r < r B 



2n + 1 



" 2^/2h B (y A r B ) n(n + 1) 



H n ( y A r B ) 



Hn (W e ) = K ( y A r B ) - ih n ( y A r B ) 

with A = L, R . The h n is a 
spherical Hankel function of 
first kind of order n, 

h{x) = x -1 /7„(x) + h;(x) , 

and R^ n , with s = e, or 
o (even or odd) are the 
spherical functions [4], [5], 

^sL(yj) = M ii n (W r ) + N sln(W r ) # x 

( 20 ) 

RsiU^r) = H^,(y R r) - N^(y R r) 

where p = 1,3, the M^' n and 
N^ n are known spherical vec- 
tor function [7]. The scat- 
tered electric field that 
comes from a LCP incident 
field or a RCP incident field 
has a similar expansion to 
(17) or to (18) 



E“(r|pJ = Q“ (r|pJ + Q“,(r|pJ = 

00 oo 

= X B n a n { L o1„(W r ) + iL< e1n(/L r )} + X { R oln (K R 0 + (W) } 



Er 8 c (r I pj = Q“, (r | P R ) + Q ^ r0 (r | p R ) + 

00 00 

= X bX^Kyj) - + X 6X{R^(//) - ^£>*0} 



Using the boundary condition (11) on r = a, we obtain [2], 
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3n ( Y L a ) h n (y R a) + j„ (y L a) h n (y R a) 
h n (y L a) h n C y R a ) + h n {y L a) h n ( y R a ) 



in (V L a ) X - in. h n (/tg) 

h n (y t a) h n (y R a) + h n {y L a) h n (y R a) 



b L in (VVO X (ygg) - -7„ (WO h„ (y R a) 
h n (V'/.a) h n {y R a) + h n {y L a) h n {y R a) 



b R in (/gg) X Mi + j„ (y R a) /;„ (yya) (2g) 

" h n (y L a) h n (y R a) + h n (y L a) h n {y R a) 

Using the asymptotic forms [4], [7], 

L'si„ (y t r) - >(n + l)(-i) n b 0 (y,r) f s L ln (r ) ( 27 ) 

RsiU/rO ~ >(n + l)(-i) n b 0 (y R r) ff ln (r) ( 28 ) 

where let us introduce LCP Beltrami angular f s L ln (r ) , and RCP 
Beltrami angular f s R ln (r) [4], satisfy by relations, 

fsi n (r) = C sln (r) + iB sln (r) , f s R ln (r) = C sln (r) - iB sln (r) (29) 

we calculate the electric far-field patterns [2], 

sc .. . . . -Z (2n + l)(-i) n_1 W n (y,r 0 ) A A ... .... .... 

g a r (r p. ) = > - — , nV L 0/ a„ {fc„(r) - if 1n (r)} (30) 

Mxr e v iw Z. 2J2n(n + 1) b 0 (y L r o ) " x elnW olnW - r v ' 



sc ^ (2n + l)(-i) n 1 H n (y R r B ) A A ... ..„ .... 

g a r (r p.) = > - — . b„ {- f p 1 n (r) - if 1 n (r) } 

^ rR/ z. 2 l 2n{n + Jj hg{ y Rrg) nt el„U 0l „K J J 
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A far-field inverse problem 

So far, all of our formu- a 7 = - i 1 + ( y L a ) 3 + 0((y L af) 

las are exact. In the asymp- I 6 

totic results to follow, L . l + fik 5 7 

there are three parameters a 2 = - i — (Y l d) + 0((y L a)) 

y A a , with A = L, R and r = a/r 0 . 

We note that the geometrical y L a^Q (34) 

parameter r must satisfy 

0<r<l because the point r ifi-R/c') 3 

source is outside of the a i = - — ttttt (V L a Y + °((Y L a T) 

sphere. I A 1 + P K ) 

We assume that |y„a| « 1, as R i(l - (3/c) 4 

well; that is we make the so- 2 36(1 + (3/c) 3 L L 

called low-frequency assump- 
tion. From (23), (24), (26), y L a^O (35) 

(27) , we obtain [2] , 



, _ 1 + [3/c a 1 2n + 1 

" 2i£ 2 (2n + 1) 



R i (1 - (3/c) n + 2 , x2n + 1 



y L a 0 



bL ~ i (1 + p/c)- 2 
" 2n£ 2 (1 - (3/c) n + 1 1 



R i(l (3/c) y2n + l 
K 2C{2n + l) {YRa) 



blL= ~ 2(l-pff ^ a ) 3+ °(^ a ) 5 ) 

b 2 = - 1(1 + P/C) 3 (y R a) 5 + 0((y L a) 7 ) 
2 36(1 - (3/c) 3 * L 



y R a -> 0 



= - i(l T (3/c) (y R a) 3 + 0 ((y R a) 5 ) 
6 

(K R a) 5 +0((y R a) 7 ) 



y R a -> 0 



Vr 9 -> 0 ( 33 ) in order to calculate g*® 

where with an error of 0((y„a) 4 ) we 

only need the following [4] 

<■„= l-3-5 - (2n - 1) = (2n)\/(2 n n\) . 

In particular, 
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5(1 1 2^ T2[f - (?) “ ] > + 0((/,a) 4 ) 



where m A = 



-1, A = L 
1, A = R 



and if A = L,R at that case A C =R,L. 



Now the scattering cross-section, by LCP or RCP spherical 
Beltrami fields, from (16), is given by the relations [2], 

< c r e = \A |Cr s (r I P jf + 4 |g~, (r I P jf] cfs(r) = 



= (na 2 ){ (1 + |3/c) [j#r 4 



[f- r 4 + (y,a) 2 (M ^ + i 6 r 6 ) + (y L a) 4 (M + ^ r 4 )] 



(1 - (3/c) 4 3 



[- r + (y L a) (3r + 



5(1 - (3/c) 2 



(1 + (3/c) 2 4 16(1 + p/c) 2 



+ (Via) 4 (3 + 



4(3 + 9(1 + pff r4) ] } + °((^a) 6 ) , / t a -> 0 



0 *%, = f A fe ( ? I Pr )\ 2 + 4 fe (r I P R )| 2 ] rfs(f) = 



= (na 2 ){ (1 6 ^ /C) r 4 + (y R a) 2 (^ r 2 + 16r 6 ) + (y R a)\^- + ^ r 4 ) ] + 

(1 + (3/c) 4 3 4 , 2 ,_ 2 5(1 + (3/c) 2 

(1 - (3/c) l 4 ' K 16(1 -(3/c) 2 ' 



+ (/ R a) 4 (3 + 



4(3 + 9(1 -$k) 2 t4) ] } + °^ a ) 6 ) ' V rS 0 



In the special case r 0 ^- oo (r -+■ 0) , by the relations (46), 
(47) we obtain [2], 



. (1 + (3/c) 2 3(1 -(3/c) 4 

3 + (1 + p/c) 2 



} (Vi a ) 4 + 0((y L a) 6 ) , y L a -> 0 
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or 

cr = (tt a 2 ) (f«a) 4 + 0((/ fi a) 6 ) , y R a -> 0 (49) 

likewise in the case y L a ->■ 0 and y R a -+■ 0 , by the relations 
(46), (47) we obtain [2], 



< re = - f A {\ 3, /c) (na 2 ) (a/r 0 ) 4 



where A = L, R , with f,((3, /c) = 



(1 - ^(3/c) 2 ; 3(1 + ^(3/c) 4 
3 (1 - m A [ 3/c) 2 



Choose a Cartesian coordi- 
nate system Oxipz, and five 
point-source locations, name- 
ly (0,0,0), (1,0,0), (0,1,0) 

( 0 , 0 , 1 ) and ( 0 , 0 , 21 ), which 
are at (unknown) distances 
r 0 , r i i r 2 > r 3 and r 4 , respec- 
tively from the sphere's cen- 
ter. The parameter 1 is a 
chosen fixed length. For each 
location, measure the leading 
order term in the low- 
frequency expansion of the 
scattering cross-section. 

Thus, our five measure- 
ments are [2], 

fflj .= ^f 4 ((3, k)na\f-Y (52) 

J = 0, 1, 2, 3, 4 

Dimensionless quantities 
related to m j are 



v = _J_ = 2 ! 1 1 (H) 2 

Vj ^ if A (P,k)TT a ( a> (5J 

j = 0, 1, 2, 3, 4 

equivalently, we obtain [2], 



r] = | y'f,((3, /c)tt /j 
j = 0, 1, 2, 3, 4 



(54) 



There are six unknowns 
namely r 0 , r lt r 2 , r 3 , r 4 and a. 
Furthermore, r 0 , r 3 and r 4 are 
related using the cosine 
rule, r 2 + r 0 2 = 2r 3 2 + 2I 2 . So, we 
can find the six unknowns. 
The center of the spherical 
scatterer is obtained from 
the intersection of the four 
spheres centers at (0,0,0) 
(1,0,0), (0,1,0) and (0,0,1), 
with corresponding radius r 0 , 
r 1 , r 2 , r 3 respectively. 
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